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Abstract 

We study the RPA equations in their most general form by taking the matrix 
elements appearing in the RPA equations as random. This yields either a unitary 
or an orthogonally invariant random-matrix model that does not appear in the 
Altland-Zirnbauer classification. The average spectrum of the model is studied with 
the help of a generalized Pastur equation. Two independent parameters govern the 
behaviour of the system: The strength a 2 of the coupling between positive- and 
negative-energy states and the distance between the origin and the centers of the 
two semicircles that describe the average spectrum for a 2 = 0, the latter measured 
in units of the equal radii of the two semicircles. With increasing a 2 , positive- and 
negative-energy states become mixed and ever more of the spectral strength of the 
positive-energy states is transferred to those at negative energy, and vice versa. The 
two semicircles are deformed and pulled toward each other. As they begin to overlap, 
the RPA equations yield non-real eigenvalues: The system becomes unstable. We 
determine analytically the critical value of the strength for the instability to occur. 
Several features of the model are illustrated numerically. 

Key words: Random-matrix theory; Random Phase Approximation 
PACS: 21.10.-k, 21.60.-n 



1 Introduction 



Random-matrix theory (RMT), introduced into physics by Wigner in the 
1950s, has found an ever increasing range of applications. In particular, it has 
been widely used to describe spectral fluctuation properties of many-body 
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systems (for a review see Ref. ppj). The main ingredients of RMT are gen- 
eral symmetry properties and the assumption that no state in Hilbert space 
plays a preferred role. The success of this approach can be ascribed to the 
fact that many physical systems are rather featureless and to a large extent 
chaotic [2]. Most theoretical work in many-body systems has, on the other 
hand, been devoted to non-statistical aspects, in particular to identifying and 
describing collective features. A prominent example is that of the giant col- 
lective states in atomic nuclei. An extreme treatment is that of the so-called 
schematic model [3] : A special (separable) form of the (residual) interaction is 
responsible for the collective features and a particular (collective) state plays 
a distinct role. A widely used many-body approximation scheme used in this 
context is the Random-Phase Approximation (RPA) first introduced to de- 
scribe the collective plasma oscillations of an electron gas. The RPA describes 
the response of the system to a time-dependent field and can be derived, for 
instance, with the help of a quasi-Boson approximation within the equation- 
of-motion method jl] or the time-dependent Hartree-Fock equations [5]. 

This paper is the first of a series of two that are written in an attempt to 
bridge the gap between the two extremes of a purely statistical (and, thus, 
"democratic") and a highly special dynamical description (the interesting fea- 
tures reside in a single state). We define a random-matrix approach to the 
RPA equations by taking the residual interaction matrix elements appearing 
in the RPA as random. This may look like an oxymoron: The RPA is designed 
to describe collective features whereas RMT is fundamentally "democratic". 
However, in the second paper of this sequel (Ref. [6], in preparation) collective 
features will be implemented in our model, and the present study of the RMT 
approach to RPA equations is a prerequisite for that work. Our paper is not 
the first to address the interplay between collective features and the back- 
ground of "chaotic states", theoretically as well as experimentally (see, for 
instance, Refs. [7] for Gamow-Teller transitions, Ref. [S] for giant resonances, 
Ref. [H] for analogue states and Ref. [TO] and the literature therein on door- 
way states) and must be seen in this context. To the best of our knowledge, 
however, this paper is the first where the problem is addressed within a purely 
random version of RPA. 

Our random-matrix model possesses the symmetry of the RPA equations. 
We study the average spectrum of the model in the limit of large matrix 
dimension with the help of a generalized Pastur equation [TTJ and compare 
the results with numerical simulations. We are especially interested in the 
instability of the RPA equations. That instability occurs when one eigenvalue 
becomes imaginary as we contiuously vary a coupling parameter. We derive 
a general criterion for the instability to develop. Changes of the coupling 
parameter cause the system to go through four different "phases" : i) the entire 
spectrum is on the real axis, ii) one part is on the real axis and the other on 
the imaginary axis, iii) one part is on the real axis, another on the imaginary 
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axis and a third one in the complex plane, and iv) the entire spectrum is on 
the imaginary axis. We are not aware of other systems, or models, showing 
this rich behaviour. 



The matrix governing the RPA equations is non-Hermitean and so is, there- 
fore, our random-matrix model. A similar model (but with restrictions that 
do not allow the eigenvalues to become imaginary) was studied in Ref. |12j . 
Our model does not belong to the ten canonical Altland-Zirnbauer ensembles 
(see Ref. p2]). We believe that, quite aside from the RPA equations, our work 
is of interest clS db CclSG study in non-Hermitean random-matrix models. 

To make the paper self-contained and to define the notation we start out with 
a brief resume of the derivation of the RPA equations in Section [2j In Sec- 
tion[3]we review the symmetries of these equations and draw conclusions about 
the distribution of eigenvalues. In Section H] the random-matrix approach to 
RPA is defined. In Section we derive the equations for the spectral density 
using the average Green's function (Pastur equations). In Sections E] and [7] 
properties of the solutions of the Pastur equations are discussed. These are 
illustrated in Section [8] with numerical simulations. In Section [9] we establish 
the critical value of the coupling strength at which instabilities occur. Sum- 
mary and conclusions are given in Section [TUJ A technical detail is deferred to 
the Appendix. 



2 RPA Equations 

For the sake of completeness, we briefly recall the derivation of the RPA equa- 
tions for particle-hole pairs of Fermions. The derivation uses the equations- 
of-motion approach. We follow Rowe's book [1]. It is essentially assumed that 
pairs of Fermions can approximately be considered as Bosons. 

We thus consider N Bosonic single-particle states with creation operators B\ 
and annihilation operators B k where k = 1, . . . , N. These fulfill the commu- 
tation relations 



We write the Hamiltonian in a form which does not conserve particle num- 
ber and which is fully analogous to the Hartree-Fock-BCS Hamiltonian for 
Fermions, 



[B k ,Bj]=S kl . 



(1) 



H = £ AlBlB, + - £ C kl B\B\ + - £ CfrBM . 



(2) 
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The matrices A and C both have dimension N. For the matrices A and C 
we consider two options. The matrix A may be Hermitean or real symmetric. 
Since the Bs commute, the matrix C is symmetric and either complex or real. 
Thus, 



A = (A°y , C = C T ("unitary case"), or 

A = (A )* = {A°) T , C = C T = C* ("orthogonal case") . (3) 

We observe that for Fermions, C in Eq. (j2j) would be antisymmetric. Let 
\RPA) denote the ground state of the Hamiltonian H in Eq. (j2J) with energy 

H\RPA) = E \RPA) . (4) 

We postulate that the excited states of the system are created by applying 
the operators 

<2l = E [X+Bl - n U B k ) (5) 



to the ground state \RPA), 

HQl\RPA) — E U \RPA) . (6) 

Multiplying Eq. (j3J with Ql and subtracting the result from Eq. (jBJ) we get 
[H,Qt]\RPA) = (E u - E )\RPA) . (7) 



The ground state is assumed to be annihilated by the operators Q v and by 
the operators 

SQ P = E (sX^B h - 5YfB\ 



obtained by arbitrary variations of the Q v , 

5Q p \RPA) = for all p . (9) 

Multiplication of Eq. (J7D from the left by SQ p and use of Eq. (jUj) yields 

[SQ p , [H,Ql]]\RPA) = (E v - E )[8Q p ,Qt\\RPA) . (10) 
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Working out the commutators, using the fact that the coefficients 8X^* and 
8Y£* are completely arbitrary, and using the first of Eqs. ([3]), we obtain 

n°x uT = (E v - E )X uT . (11) 



We have defined 

H ° = ( A n* _ Fag\* ) ( 12 ) 



-C* -(A )* 
and 

X u = {Xl,Y»). (13) 



Eqs. (llip to ({TBI are the RPA equations for Fermionic particle-hole pairs. 
Except for one important difference these equations look very similar to the 
Hartree-Fock-BCS equations for Fermions: The matrix C is symmetric while 
in the BCS case, the pairing potential is antisymmetric. As a consequence the 
BCS Hamiltonian is Hermitean while the matrix 7i° in Eq. ffl2|) is not. This is 
why the 2N x 2N matrix 7i° does not belong to one of the symmetry classes 
studied in Ref. [T3l. 



3 Spectral Properties 



Since the RPA matrix 7i° is not Hermitean, the eigenvalues need not be real. 
We investigate some of the consequences. 

The matrix 7i° in Eq. (fT2l has the following two symmetry properties. With 
1 7v the unit matrix in N dimensions and 



Mi 



1 N 
In 



m, = (\ n .J..) , (in 



N 



we have 



M^YM^-H , 

M 2 (nyM 2 =n° . (is) 
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Let be an eigenvalue and ^ be an eigenfunction of the RPA matrix, 
H°V, = E^n . (16) 

From the first of Eqs. (IT5]) and (Mi) 2 = 1 2 n it follows that (-E*) is also 
an eigenvalue with eigenfunction A4i$f* From the second of Eqs. (IT5|) and 
(M.2) 2 = I2JV it follows that E* is also an eigenvalue with left eigenfunction 
We conclude that with also E* —E^ and — E* are eigenvalues of the 
RPA matrix US ' . In other words, real and purely imaginary eigenvalues come 
in pairs with opposite signs, while fully complex eigenvalues with 3?_E M 7^ 
and SsE^ 7^ come in quartets that are symmetric with respect to both the 
real and the imaginay energy axis. 

For C = the matrix 7i° is Hermitean or real symmetric. Then the eigenvalues 
of 7i° are real and appear in pairs with opposite signs. While level repulsion is 
the hallmark of a Hermitean Hamiltonian, we show that the non-Hermitean 
part C of H° leads to repulsion (attraction) of levels with the same (oppo- 
site) signs, respectively. With increasing strength of C level attraction leads 
to coalescence of pairs of eigenvalues with opposite signs at energy E — 0. Co- 
alescence is a signal for the instability of the RPA approach. With a further 
increase of the strength of C the two coalesced eigenvalues leave the point 
E = in opposite directions along the imaginary E-axis. On the imaginary 
axis, a pair of eigenvalues may coalesce again and then move away from the 
imaginary axis in opposite directions along a line that is parallel to the real 
.E-axis. Because of the symmetry of the RPA matrix, such behavior must 
occur in parallel on the positive and on the negative imaginary axis. 

Before giving the general argument, we demonstrate these statements using 
two simple examples. To be specific we consider the unitary case. First, we take 
N = 1 and label the matrix elements in an obvious way by a and c, respectively. 
Here a is real and c is complex. The eigenvalue equation is quadratic with the 
two solutions Q = ±-y/a 2 — |c| 2 . Increasing the value of |c|, we pull the two 
eigenvalues toward each other. The two eigenvalues coalesce at E — for 
|a| = |c| and become purely imaginary for |a| < |c|. Second, we take N = 2. 
We assume that A has been diagonalized with eigenvalues E\ and E 2 where, 
in the spirit of the RPA, we assume < E\ < E^. The complex matrix 
elements of C are labelled ci,C2,cu in an obvious way. The secular equation 
is 

(n 2 - E 2 )(n 2 - e\) + \c x \ 2 ($i 2 - El) + \ C2 \ 2 (n 2 - eD 

+2|c 12 | 2 (fi 2 - E X E % ) + I det(C)| 2 = (17) 
which we write in the form 
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ft 4 - n 2 (T A - T c ) +D C + E\E\ 

= El\c 2 \ 2 + E 2 2 \ Cl \ 2 + 2E 1 E 2 \c 12 \ 2 

= T AC . (18) 

Here T A = E\ + E 2 , T c = Trace[CC*], D c = \ detCf , and the last part of 
Eq. (jTSJ) defines T AC . This quadratic equation for Q 2 has the solutions 

1 



n 2 = -{T A - T c ) ± -y/(T A - T c ) 2 + 4T AC - 4E 2 E 2 - AD C . (19) 



The argument of the square root vanishes for 

D c = T AC + \{T A - T c ) 2 - E\E\ (20) 

while two eigenvalues coincide at E = if 

D c = Tag ■ (21) 

By definition we have T A > 0, Tc > 0, T^c > and Dc > 0. Inspection 
shows that also (T A — Tc) 2 — 4Ef _Ef > 0. As we turn on the coupling matrix 
C, Tc,T A c, and Dc grow monotonically from zero, and the condition (j2ip 
is met prior to condition (120]) . This shows that the smaller of the two posi- 
tive eigenvalues coalesces with its negative counterpart at zero (and becomes 
imaginary) without the two positive (or the two negative) eigenvalues ever 
coalescing. 

We show that coalescence of two eigenvalues of opposite signs at E = (with- 
out coalescence of pairs of positive or pairs of negative eigenvalues) is a general 
property of the RPA equations. We denote the orthogonal projectors onto the 
two iV-dimensional subspaces appearing explicitly on the right-hand side of 
Eq. (fl2i) by Qi and Q 2 , respectively. With E the energy, we eliminate the sub- 
space with projector Q 2 at the expense of introducing into the first subspace 
the effective Hamiltonian H = A -C{E + A *)' 1 ^ . The term added to A is 
a Hermitean (or a real symmetric) operator and, thus, causes level repulsion 
between the eigenvalues of A which, therefore, cannot coalesce as the matrix 
elements of C increase in magnitude. (The non-Hermiticity of the matrix 7i° 
leads to the negative sign of the term —C(E + A°*)~ 1 C*, while a Hermitean 
H° would have given a positive sign. But level repulsion is independent of the 
signs of the matrix elements connecting the states). Going to the eigenvalue 
representation of the matrix A and assuming that all eigenvalues are positive, 
we find that the trace of —C(E + A *)~ 1 C* is negative for all E > and 
largest for E = 0. This shows that with increasing strength of C the center of 
the spectrum of the positive eigenvalues is shifted towards smaller values, the 
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shift being strongest for the levels near E = 0. We conclude that the interac- 
tion described by C causes level attraction between the eigenvalues of A and 
those of —A * while it causes level repulsion amongst the eigenvalues of A 
and amongst those of A *. 

To investigate the coalescence of eigenvalues at E = we follow Kato's 
book [T5]. We denote by TC°(z) the matrix obtained from 7i° by the replace- 
ment C — » zC. We assume that the eigenvalues of A are all positive. We 
consider the characteristic polynomial P(Q,z) = det[H°(z) — QI2N] for com- 
plex values of z. The following statements follow directly from Chapter II. 1 of 
Ref. [15]: The roots Q(z) of the equation P(fl, z) = are branches of analytic 
functions of z with only algebraic singularities. At the point z = z c where two 
or more eigenvalues coalesce two eigenvalues assume the common value Q c if 
the two equations 



are both fulfilled. For three or more eigenvalues to coalesce, further equations 
involving higher derivatives of P must also be fulfilled. Such equations impose 
constraints on HP. Therefore, the generic case is the coalescence of two eigen- 
values. (This is a generalization [16] of the Wigner-von Neumann non-crossing 
rule for eigenvalues). We confine ourselves to that case. Two eigenvalues co- 
alescing at z c generically possess a square-root branch point at z c (the case 
where the two coalescing eigenvalues are holomorphic at z c also leads to con- 
straints on TC°, is non-generic and, thus, likewise excluded [16J). The same 
argument applies to the coalescence of eigenvalues on the positive and on the 
negative imaginary i?-axis. 

The behavior of two coalescing eigenvalues was already studied in the simplest 
case N = 1 introduced above. Both eigenvalues are real and have opposite 
signs for z < z c = \a\/\c\ and are purely imaginary and have opposite signs for 
z > z c . With increasing real z < z c the two eigenvalues move along the real 
axis towards each other until they coalesce at z — z c . As z increases further, 
the two eigenvalues leave the point zero and move along the imaginary axis in 
opposite directions. 

That behavior is not confined to the case N = 1 but represents the generic 
situation. We consider a point of coalescence z c with z c real where the two 
coalescing eigenvalues (denoted by Oi and f^) both have the value zero. We 
assume that for z real and z < z c , Q 1 (fi 2 ) is positive (negative, respectively). 





and 




(23) 
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For z in the vicinity of z c we therefore must have 



(24) 



where (3 is some positive constant. Eq. (1241) shows that while f2i (f2 2 ) is positive 
(negative) for z real and z < z c , both Q± and are purely imaginary and 
carry opposite signs for z real and z > z c . 

The case of two eigenvalues coalescing on the imaginary axis is similar. With 
z c the critical strength for coalescence at Q = Q c with Q c purely imaginary, 
the two eigenvalues are purely imaginary for z < z c . Therefore, Eq. (I24p takes 
the form S7x,2 = ± (3y/z c — z where j3 is purely imaginary. It follows that 
for z > z c the two eigenvalues leave the imaginary axis in opposite directions 
along a straight line that runs parallel to the real E-axis. 



4 Random— Matrix Approach 



In order to study the generic properties of the RPA equations, we consider 
a Gaussian random-matrix model. In Fermionic systems, the RPA equa- 
tions (fTTj) to ffl~3l) account for the particle-hole interaction. As that interaction 
tends to zero, the matrix C vanishes while the matrix A tends to a diagonal 
matrix with elements given by the single-particle energies of the particle-hole 
pairs. We assume that all these energies have the same value which we denote 
by r. In the nuclear case we have r = ftw where u is the frequency of the 
harmonic-oscillator potential. We write A in the form 

A = rl N + A . (25) 



The Hermitean (or real symmetric) N-dimensional matrix A and the complex 
(or real) symmetric matrix C represent the particle-hole interaction. We as- 
sume A to be a member of the Gaussian unitary ensemble (GUE) or of the 
Gaussian orthogonal ensemble (GOE), as the case may be. The matrix ele- 
ments are complex or real Gaussian-distributed random variables with mean 
values zero and second moments given by 

A 2 

(A^ u A pcr ) = —S^5 up (unitary case), or 
A 2 

(A^Apu) = —(5 pa 5 up + 5 w 5 ua ) (orthogonal case) . (26) 

The indices run from 1 to N . The angular brackets denote the ensemble av- 
erage. For iV ^> 1, the average GUE and GOE level densities have the shape 
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of a semicircle with radius 2A. The GUE (GOE) is invariant under unitary 
(orthogonal) transformations A — > UAU^ (A — > OAO T ) where UW = In 
(OO t = In, respectively). The elements of C are assumed to be Gaussian- 
distributed random variables with zero mean values and second moments 



(C^ u C* a ) = jj(5^Jvp + S w S ua ^ , {C^ v C pa ) = (unitary case) , 

(C^Cpa) = (6^6^ + 5p P 5ua^j (orthogonal case) . (27) 

The ensemble of matrices C is invariant under the transformations C — > 
UCU T where UW = In and where T denotes the transpose. The elements of 
A and of C are assumed to be uncorrelated. The spectrum of 7i is expected to 
be real as long as the matrix elements of C are not too large. This is why we 
parametrized the Gaussian distribution of C in terms of an extra parameter 
7. We simplify the presentation by considering only the unitary case in the 
sequel. The orthogonal case is obtained with a slight change of notation; our 
conclusions apply to both cases. 

It is useful to display the statistical properties of A and of C in some detail. 
Writing the elements of the matrix A in the form 

V = 3LV + *3 V (28) 



we note that dtA^ is real symmetric and QA^ is real antisymmetric. The 
elements of $t,Ap U and of are uncorrelated Gaussian-distributed random 
variables with zero mean value and second moments given by 



(29) 

Similarly, writing the elements of the matrix C as 

CV = KC V + ^C M „ (30) 



we note that both 3?C M ^ and ^sC^ are real symmetric. The elements of ^tC^ 

correlated Gaussian-distri 
zero mean value and second moments given by 



and of ^sC^ are uncorrelated Gaussian-distributed random variables with 



<(^) 2 > = (1 + V)^, 
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((^) 2 > = (1 + <Wt^- (31) 
With these definitions, the RPA matrix 7i° takes the form 



rl N \ / A C 

-rl N ) + \-C* -A* 

rl N 

-rl N 



+ H . (32) 



We refer to the joint invariance properties of A and C as to "generalized 
unitary invariance" ("generalized orthogonal invariance", respectively). The 
generalized unitary invariance says that the ensemble of matrices TC is invariant 
under the transformation 



It was mentioned already that since TC is not Hermitean, the random-matrix 
ensemble constructed above does not belong to one of the classes listed in 
Ref . [?] . Non-Hermitean or Non-Cartan random-matrix models have recently 
received some attention, see Refs. [T4|fT7] . It seems that our model appears in 
Table 7 of Ref. [14J as class 21b (lower sign). We are not aware, however, of 



any detailed analysis of the model. 

We study our random-matrix models in the limit iV — > oo. In that limit, the 
parameter N disappears, and each of our random-matrix models is character- 
ized by two parameters. We introduce the first by noting that for 7 = the 
average spectrum consists of two semicircles. Each of these has radius 2A. We 
measure all energies in units of that radius. The first parameter r is then the 
distance of the center of each semicircle from the origin at E — 0. The second 
parameter measures the strength 7 of the matrix C in units of the strength A 
of the matrix A, see Eqs. (|26l) and (1271) . It is given by 



« 2 = £ ■ (34) 



We will use these parameters when we treat the Pastur equation. Using the 
results of Section [3] we anticipate that with increasing a 2 , the two branches 
of the average spectrum are deformed and move toward each other. This is 
indicated schematically in Fig. [TJ 

The average spectrum of our random-matrix model (|32|) takes a simple form 
in two limiting cases. Aside from the simplification encountered for C = 
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Fig. 1. Expected deformations and shifts of the two branches of the average spectrum 
with increasing a 2 (schematic). 

and discussed above, another simplification arises in the opposite extreme 
case when the spectrum is dominated by the matrix C. To see this we put 
in Eq. (1521) r = 0, A = and C = iD and write the eigenvalues as iQ. The 
resulting eigenvalue equation is 

U°. o)* = n *- (35) 



That equation possesses chiral symmetry. It is straightforward to derive the 
Pastur equation and to see that the average spectrum has the shape of a 
semicircle centered at zero with radius 2A. We conclude that as a 2 in the RPA 
matrix in Eq. (|32|) increases from zero to very large values, the eigenvalues 
move from the two semicircles on the real axis onto a single semicircle on the 
imaginary axis. As they do so, they may intermittently leave both the real 
and the imaginary i?-axis and form groups of quartets arranged symmetri- 
cally with respect to both these axes. This is exemplified in Fig. [2J For r = 4 
and several values of a 2 , we show the distribution of eigenvalues in the com- 
plex energy plane. These were obtained by diagonalization of randomly drawn 
matrices of dimension 2N = 40. We see that before all eigenvalues reach the 
imaginary axis at a 2 = 150, a quartet is formed at a 2 = 50. We have found 
numerous such quartets in other random realizations. Our simulations also 
suggest that in the unitary case, the semicircle formed for very large values 
of a 2 on the imaginary axis has a gap at E = 0. In certain chiral ensembles, 
the average level density is known [i~8][T9] to vanish at E = 0, and our result 
probably relates to this fact. We have not followed this up any further. 
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a=50 



cf=150 




Fig. 2. Evolution of the location in the complex energy plane of the eigenvalues for 
r = 4 and 2N = 40 when the coupling is increased, i.e. for a 2 = 1, 30, 50, 150. 

5 Pastur Equations 



We ask: How do the matrix elements of C affect the average RPA spectrum? 
Which is the smallest value of 7 where the average RPA spectrum becomes un- 
stable? We answer these questions by calculating the average retarded Green's 
function of the system. It is given by 

( G ( B )) = ((£+1 2 „-C? -rl J -")"'}■ (36) 

For sufficiently small values of 7 2 , we expect the average spectrum of H° to 
be real and to consist of two branches, one at positive and one at negative 
energies. The instability of the RPA occurs when the two branches begin to 
touch. 

To calculate the ensemble average, we expand G(E) in powers of H. This 
yields 

00 

G(E) = G (E) + £ G (E)(HG (E)) n . (37) 

Here G (E) is the Green's function of the system for H = 0, 

Go(B )=( B+ l 2w -( rl " _£))". (38, 
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For N ^> 1 the calculation of the ensemble average of G is standard: We 
average each term in the sum in Eq. (1371) separately using Wick contraction. 
For N ^> 1, we keep in each term only the contributions with the maximum 
number of free summations over the N level indices. The resulting series can 
be resummed and yields the Pastur equation 

(G{E)) = G Q {E) + G Q (E)(H(G(E))H)(G(E)) . (39) 



Using the same notation as in Section [3j we introduce the orthogonal projec- 
tion operators Qi,i = 1,2 onto the two ^-dimensional subspaces displayed ex- 
plicitly on the right-hand side of Eq. (133]) . with Q\ = Q 1 , Q\Q 2 = 0, Q1 + Q2 = 
1 2 n etc. The average Green's function is written as 

(G{E)) = £ Qi(G(E))Qj . (40) 

The terms with % 7^ j on the right-hand side of Eq. (I40p vanish. Indeed, it 
takes an odd number of matrices C and an even number of matrices C* or 
an even number of matrices C and an odd number of matrices C* to connect 
subspaces 1 and 2. The elements of C being independent Gaussian random 
variables with mean values zero, such terms vanish upon averaging. This leaves 
us with 

(G{E)) = Q 1 (G(E))Q 1 + Q 2 (G(E))Q 2 . (41) 



We project the Pastur equation (139]) onto each of the two subspaces, use 
Eq. (jHJ), observe that Gq(E) commutes with Q\ and Q2, and obtain 



Qi(G(E))Q 1 = Q 1 G (E)Q 1 + (Q 1 G (E)Q 1 ) 

x{((QinQ 1 ){Q 1 (G(E))Q 1 )(Q 1 HQ 1 ))(Q 1 (G(E))Q 1 ) 
+((Q 1 HQ 2 )(Q 2 (G(E))Q 2 )(Q 2 HQ 1 ))(Q 1 (G(E))Q 1 )} 
= QiG (E)Q 1 + (Q 1 G (E)Q 1 ) 

x{(A(Q 1 (G(E))Q 1 )A)(Q 1 (G(E))Q 1 ) 
-(C(Q 2 (G(E))Q 2 )C*)(Q 1 (G(E))Q l )} (42) 

and 



Q 2 (G(E))Q 2 = Q 2 G (E)Q 2 + (Q 2 G (E)Q 2 ) 

x{((Q2HQ 2 )(Q 2 (G(E))Q 2 )(Q 2 HQ 2 ))(Q 2 (G(E))Q 2 ) 
+ {{Q2HQ l ){Q 1 {G{E))Q 1 ){Q 1 HQ 2 )){Q 2 (G{E))Q 2 )} 
= Q 2 G (E)Q 2 + (Q 2 G (E)Q 2 ) 
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x{(A*(Q 2 (G(E))Q 2 )A*)(Q 2 (G(E))Q 2 ) 
-(C*(Q 1 (G(E))Q l )C}(Q 2 (G(E)}Q 2 )} . (43) 

We perform the averages over the matrices A and C using Eqs. (126]) and (127)) . 
define for i = 1,2 the functions 

= A T raceg,,(G( J E;))g j , (44) 



use that 

Q 1 G Q (E)Q 1 = -t^—In , Q 2 Go(£)Q 2 = -J—^, (45) 
E — r E + r 



take the traces of Eqs. (142)) and (H3l and obtain, using iV 1 



^ A w T 2 . . 

0i = p + p A(«7! - -^(x 2 )Aa 1 , 

E — r E — r X z 

XX 7 2 

a2 = eT7 + eT7 a((T2 " v ai)X(T2 ■ (46) 



With a 2 = 7 2 /A 2 these equations can be rewritten in the form 
A 



E — r — A(a"i — a 2 <72) 



^ = cu vr • ( 47 ) 

E + r — X(cr 2 — a. (J\) 

These are two coupled quadratic equations in the unknown functions <J\{E) 
and a 2 (E). We refer to both Eq. (f39l) and to Eqs. (147)) as to the (generalized) 
Pastur equation(s). A physical interpretation of the functions a%{Ei) with i = 
1, 2 is given at the end of Section [6j 



6 Properties of (G(E)) 



Before averaging, Trace G(E) can be expressed in terms of the eigenvalues of 
TC°. We confine ourselves to sufficiently small values of |^y| so that for N — > oo, 
the average spectrum of H° consists of two separate branches, the positive 
branch with eigenvalues E^ > 0, and the negative branch with eigenvalues 
— E^ < 0. Here, \i — 1, . . . , N. We have used the symmetries of the spectrum 
derived in Section [21 
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We distinguish the advanced and the retarded Green's functions by appropri- 
ate indices. With = E ±i5 where E is real and 5 positive infinitesimal, we 
have 



Trace G(E) retMv = ^ ^ 



2 N j 



, 2 N 



i=l n=l 

Eq. (T48l) is standard for Hermitean operators. But Ti is not Hermitean. There- 
fore, we derive Eq. (j4"51) in the Appendix. Averaging over the ensemble converts 
the sum over delta functions into the average level density p{E), 

, 2 N . 

EE^ + H%) - (49) 
\ i=l M =l ' 



where is normalized such that the integral over all energies yields 2N. It 
is obvious that 

p(E) = p{-E) . (50) 



The average level density is a symmetric function of E. It differs from zero 
within two energy intervals, one extending over positive energies with end 
points, say, E a and Eb where Ej, > E a > and the other, extending over 
negative energies between — Eb and —E a . 

We use p(E) to write (G(E)) in the form 

(Trace G(E) ietMv ) = J dE' ^T^7 ■ ( 51 ) 



Eq. floTl) shows that the analytic function (TraceG(-E)) has two branch cuts. 
The branch cut at positive energy extends from the branch point at E = E a to 
the branch point at E — Eb- The branch cut at negative energies is the mirror 
image of the first with respect to the origin. The discontinuity of (TraceG(-E')) 
across the cut is given by 

(TraceG(£) adv ) - (TraceG(£) ret ) = 2mp(E) . (52) 



For Hermitean problems, the retarded and the advanced Green's functions are 
related by complex conjugation. However, TL is not Hermitean, and (G(E) Tet ) 
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and {G(E) a dv) are not complex conjugates of each other. Inspection of Eqs. (|36|) 
and (132]) shows that with Ai defined in Eq. (fl4l) we have 

- M(G^ xet {-E))M = (Gadv,ret(S)) - (53) 



As long as the two branches of the spectrum do not touch, the functions 
<Ji{E) and &2{E) have a simple physical interpretation. That interpretation 
uses Eq. (jUJ) and the spectral decomposition of G(E) given in Eq. fPfBl) . We 
first use the results in Section [3] to conclude that if ^ is a right eigenfunction 
of 7i° with real eigenvalue E^ > 0, then \l/J t M 2 is the left eigenfunction to the 
same eigenvalue, while M.^*^ and ^.MiA^ are the right and left eigenfunc- 
tions, respectively, to eigenvalue — E^. We assume that the eigenvalues are 
not degenerate. We expand G(E) itself (rather than its trace) in the form of 
Eq. (l4"8j) and project the result onto the subspaces with index 1 or 2. We use 
M 2 Qi = Qi and M2Q2 = -Q2- This yields 



N 
N 

-nJ2(^l\M 1 Q 1 M 1 \%)5(E + E,) , 

N 

sg 2 G(E) ret g 1 = +7T J2i*t W^ME - 

N 

+7r^K\M 1 Q 2 M 1 \%)6(E + E fl ) . (54) 
fj,=i 

Eqs. ( |54l show that ai(E) (a 2 (E)) measures the total average spectroscopic 
strength of all states in the subspace with index 1 (with index 2, respectively). 
For a 2 = 0, 0\ [o<i) receives non-vanishing contributions from the positive- 
energy states (the negative-energy states) only. As a 2 increases, that situation 
is expected to change. This is confirmed by our numerical calculations. 

The analytic properties of (G(E)) used in the present Section may fail to hold 
when a finite fraction of eigenvalues (i.e., a fraction of order N) is neither 
purely real nor purely imaginary. In that case use of the Pastur equations as 
in the present paper may not give correct results because the Green's function 
has other singularities besides branch points, see f.i. Ref. [20] . 
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7 Properties of the Solutions of Eqs. ( 14 7j) 



We are going to use the solutions tT 12 of Eqs. (|4"T|) to determine the average 
level density p(E) of the system. By elimination of cr 2 , Eqs. ( |4"T|) can be re- 
duced to an equation of fourth order for <7i, see Eq. fl60l) below. Thus, Eqs. (jUJ) 
possess four pairs of solutions (a[ l \ ) z = 1, . . . , 4. It is neccessary to deter- 
mine which of these is relevant to our problem and how that solution relates 
to p(E). We use perturbation theory for small values of a 2 and continuity to 
establish that connection. In order to deal with non-overlapping spectra we 
take 2A < r. 

We first take a 2 = 0. Then Eqs. f|4"Tj) decay into two uncoupled equations for 
the unperturbed spectroscopic strenghth S 12 . With = (E + (— )V)/(2A) the 
equations read 

s» = ^rs; ' i = 1 ' 2 - ( 55 ) 

The solutions are 

E< = ei ± ^e? - 1 , i = 1, 2 . (56) 



The solutions have square-root singularities (branch points) at 6{ = ±1. For 
Si (S 2 ), the two branch points are at E a = r — 2A and at Eb = r + 2A (at —Eb 
and at — E a , respectively). The branch cut extends along the real £"-axis from 
E a to E b (from —E b to —E a , respectively) and connects the two sheets of the 
Riemann surface. We define the physical sheet as the one where the imaginary 
part of Ej is positive when we approach the branch cut from below the real 
axis. Then the square root carries a negative sign for E real and E ^> r, and 
Ej ~ X/E for E real and \E\ ^> r. On the second sheet we have Ej rs E/X 
for E real and \E\ 3> r. For both solutions, the level density differs from zero 
only if E lies on the branch cut, see Eq. f)52p . In both cases, the average level 
density is given by the value of 3Ej on the physical sheet and has the form of 
Wigner's semicircle law (N/(n\))(l - [(E - (±)r)/(2A)] 2 ) 1 / 2 . 

We now consider the case of small a 2 and solve Eqs. (|47|) perturbatively up 
to first order in a 2 . This yields 



a 1 = E 1 -« 2 [E 1 ] 2 E 2 , 

a 2 = E 2 -a 2 [E 2 ] 2 E 1 . (57) 

We note that the end points of the spectrum are not changed. We would 
expect that with increasing a 2 , the branch points are also shifted. Such a shift 
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is, however, beyond the reach of a straightforward perturbative approach, see 
the end of this Section. We focus attention on o\ and consider the physical 
sheet of Ei where 3Ei > for E a < E < E^. Because of the symmetry 
of Eqs. (I5TI) our arguments carry over straightforwardly to cr 2 . The coupling 
between states in subspaces 1 and 2 shifts some of the spectroscopic strength of 
subspace 1 into subspace 2 and, thus, into the interval —Ef, < E < —E a , and 
vice versa. Both this gain and the corresponding loss of strength in the first 
interval are given by the last term in the first of Eqs. (I5T1) . The first of Eqs. (15^1) 
shows that (whatever choice we make for 3<7i) the contributions to 3cri from 
the positive- and from the negative-energy states must have the same signs. 
We have chosen 3Ei > for E a < E < Ef,. Therefore, in the interval — Ef, < 
E < —E a , the imaginary part of o"i given by — S([Ei] 2 E 2 ) = — [Ei] 2 SE 2 
must be likewise positive. That means that 3E 2 must have a negative sign, 
or that we must choose E 2 to lie on the second sheet. We conclude that Sa' 2 
must be negative: The physically meaningful solutions of Eqs. ([57]] are the 
ones where ^sa± and 3<t 2 carry opposite signs. In the interval E a < E < Ef,, 
the imaginary part of the last term in the first of Eqs. floTl) has the form 
— 2a 2 9ftE 1 9E 1 E 2 . Here E 2 « E/X is positive on the second sheet and grows 
monotonically with E, the positive function 3*Ei is given by the semicirle law 
and is symmetric about E = r, and KEx is antisymmetric about that point 
and negative for E < r. It follows that — 2a 2 9?EiQEiE 2 is positive (negative) 
for E < r (E > r, respectively). Moreover, for equal distances of E from r the 
negative part has larger magnitude than the positive one. Aside from its loss 
to the negative energy interval, the level density in the positive energy interval 
is shifted toward smaller energies. This was expected, see Section [3J The total 
balance of all these changes must be zero, and the energy integrals over Qa± 
and QEi must be equal since / p(E) = 2N by definition as long as the two 
branches of the spectrum do not touch. We have not checked that statement 
analytically for the terms in Eqs. flBTl) but have tested our numerical results 
accordingly, see Section [SI 

In summary we have shown that for small a 2 we must choose those solutions 
of Eqs. (jUJ) for which the imaginary parts of 0\ and <r 2 have opposite signs. 
Using continuity we conclude that this statement applies for all values of a 2 . 
Eqs. (jH]) being real, the non-real solutions (<7i, cr 2 ) come in complex conjugate 
pairs. Without loss of generality we can choose the solution for which Qa± is 
positive. For E > the total level density of the system is then given by 

N 

p(E) = +—%a 1 (E) + a 2 (E)) . (58) 



We observe that Eq. (!58l) and the symmetry of the system imply that p(E) 
vanishes at E = 0. This is true even when some of the eigenvalues have moved 
off the real axis. 
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To investigate the motion of the branch points due to changing a 2 , we use a 
modified version of perturbation theory. We substitute in the first of Eqs. (1471) 
for o"2 the asymptotic value E/X on the second sheet and obtain 

ai = WTTZ 2^ • ( 59 ) 

E(l + cr) — r — o\ 



Solving this quadratic equation we find for the new branch points E' a and E' h 
the values E' a = E a /(1 + a 2 ) and E' b = Eb/(1 + a 2 ). That shows that both 
branch points are shifted toward smaller energies as expected. 



8 Numerical Results 

By elimination of 02, we obtain from Eq. (HTj) an equation for <j\. It reads 
aa\ + ba\ + co\ + da x + e = (60) 



where 



a = A 2 (l-« 4 ) , 

b = X{E-r){a 4 -2)-Xa 2 {E + r) , 
c= (E-r) 2 + a 2 (E 2 -r 2 ) + 2X 2 , 
d=-2X(E - r) - Xa 2 (E + r) , 

e = A 2 . (61) 

Without loss of generality we put 2A = 1. Equivalently, we go over to dimen- 
sionless variables by measuring all energies in units of 2A. Then the semicircles 
have equal radii of value unity, r is, for a 2 = 0, the distance between the center 
of each of the two semicircles and the origin, a 2 the strength of the coupling, 
and E the dimensionless energy. 

Except for the non-generic cases where the coefficients a and b are equal to 
zero, Eq. (160]) has four solutions denoted by (i = 1, 2, 3, 4). We may either 
have 4 complex solutions (all four imaginary parts do not vanish identically), 
or 2 complex and 2 real solutions, or 4 real solutions. Complex solutions always 
come in complex conjugate pairs. To all four solutions corresponds a unique 
value of o"2 obtained from the first (linear) equation of Eqs. (T4T|) . If <j\ is 
real (complex) then 02 is also real (complex). Thus we have the same three 
possibilities for the pair (er^ , a^) '■ 4 real, 4 complex or 2 real and 2 complex 
pairs. As explained in Section [7J we look for the complex solution with S(eri) > 
and 3(cr 2 ) < 0. It turns out that that solution occurs for the cases when we 
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have 2 real and 2 complex pairs of solutions, and when we have four complex 
solutions. 



We summarize the properties of the solutions of Eqs. (14"T|) found numerically 
in Fig. [3j For three values of r, we show in the (E, a 2 ) plane four domains 
labelled (I) to (IV). These symbols stand for 

• (I) 4 complex pairs, 

• (II) 2 real and 2 complex pairs with equal signs for £y(o"i) and ^(c^), 

• (III) 2 real and 2 complex pairs with opposite signs for S(o"i) and ^(c^), 

• (IV) 4 real pairs. 

The physically interesting solutions lie in domains I and III. In all three cases 
shown in Fig. [3l the border lines between domains III and IV are not verti- 
cal but slightly inclined: With increasing a 2 , the spectra are pushed toward 
E = 0. Similarly, the domains II and IV located at the center become squeezed 
and finally disappear. For a 2 > 1 we only have 4 real solutions (case IV) or 2 
real and 2 complex solutions whose imaginary parts have opposite signs (case 
III). 





(III) ~ 




(IV) 


/ (I) \ 


(IV) 










/ H \ 






/ '* " 1 \ 









(IV) 










/ < n » \ 








(III) 




(in) 




(IV) 








(IV) 




/ (I) 




(i) \ 








•' (IV) \ 







(a) 



-1 o 
E 



(b) 



1 






(IV) 










(III) 




(m) 






(IV) 








(IV)- 


B 

0.5 


















/ / '•. (II) \ 










/ J) 


(IV) 


d)\ 








/ 









(c) 



Fig. 3. (Color online) The four domains (I) to (IV) in the (E, a 2 ) plane for three 
values of r. (a): r = 1.2, (b): r = 2.0, (c): r = 4.0 (See text for further explanation). 

In Fig. H] we compare for E > the average level density obtained from the 
Pastur equation with the actual level density obtained by numerical simula- 
tion. We have constructed 100 realizations of the RPA matrix in Eq. (T2"Tj) with 
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dimension 2N = 100 for the unitary case for r = 2 and three values of a 2 by 
drawing the elements of the matrices A and C from Gaussian distributions. 
We have determined the eigenvalues by diagonalization. The results are shown 
as histograms. The imaginary parts of o\ and cr 2 are given as dot-dashes and 
dotted lines, respectively, their sum as a solid line. The very close agreement 
shows that the Pastur equation is capable of predicting average spectra quite 
accurately even for matrices of fairly small dimension. 

r=2, <x 2 =o 



0.5 



- - Im(o,) 



1 1.5 
E 



2.5 



r=2, a =1 




r=2, a =4 




Fig. 4. (Color online) The average level density p(E) versus E for E > as obtained 
from the Pastur equation. Black dashed-dotted line: Qui, red dashed line: 9cJ2, blue 
solid line: sum of the two. For comparison, histograms corresponding to a numerical 
simulation with 100 RPA matrices H of dimension 2iV = 100, r = 2 and different 
values of a 2 . 
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Fig. 5. (Color online) The average level density p(E) versus E for E > as obtained 
from Eq. (|58|) . Upper (lower) panels: p(E) for r = 1.2 (r = 4.0) and for several values 
of the coupling parameter a 2 . Dahed-dotted, dashed and full lines as in Fig. [H 
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Fig. [5J shows the level density p(E) as obtained from Eq. (1351) for E > 0. 
For a 2 = we obtain two semicircles. These are separate for r > 1 but 
would overlap for r < 1. With increasing a 2 , the semicircles become ever 
more deformed and are pushed towards each other, as expected from the 
perturbative result in Section [71 Average spectra that are separate for r > 1 
and a 2 = eventually touch. We observe that the critical value of a 2 where 
this happens (denoted by a 2 rit ) increases with increasing r, starting out from 
a crit = at r = 1. 

In Fig. [6] we show the value of the integral I = (1/ (2iV)) / dE p(E) versus a 2 . 
For a 2 = the integral has the value unity. For < a 2 < a 2 Tit the value of 
/ should be very close to unity. In that regime, some elements of the matrix 
C for some rare realizations may have large values so that some eigenvalues 
become non-real and do not contribute to I. The resulting decrease of I is 
expected to vanish as iV — > oo, however, and is not visible in Fig. [61 For 
the value of I decreases monotonically with a 2 because an ever 
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increasing macroscopic fraction of eigenvalues leaves the real E-&xis. 
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Fig. 6. The norm integral (l/(2N))p(E) f dE versus a 2 . See text for further expla- 
nation. 



9 Critical Strength 



The critical strength a 2 Tit is defined as the smallest value of a 2 for which for 
two solutions whose imaginary parts carry different signs, So"i and 3o"2 do 
not vanish at E — 0. In order to find an analytical expression for a 2 rit as a 
function of r, we put in Eq. ( |60l) E = and define, using Eqs. ( IBTl) . the new 
variable y by 



(62) 
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From Eq. (!60|) we obtain 
y 4 + uy 2 + v y + w = 



(63) 



where 



3b 2 c 




8a 3 2a 2 a ' 



/ b \4 c / b \ 



bd e 
4a 2 a 



(64) 



Eq. (163]) is an equation of fourth order. The solutions are given by a combi- 
nation of the three roots of the third-order polynomial R(Z), 



The polynomial R(Z) has real coefficients. Therefore, the following possibilities 
exist. 

• (i) The three roots of R(Z) are real and have equal signs; the quartic equa- 
tion fl63l has 4 real solutions. 

• (ii) The three roots of R(Z) are real, only two roots have the same signs; 
the quartic equation ([6"3"j) has 4 complex solutions. 

• (iii) Two roots of R(Z) are complex and one is real; the quartic equation f[6"31 
has 2 real and 2 complex solutions. 

The roots of R(Z) are determined via the Cardan method. The substitution 



R(Z) = Z 3 + 2uZ 2 + (u 2 - Aw)Z - v 2 . 



(65) 




(66) 



eliminates the quadratic term in Eq. fl65|) and yields the equation 



X 3 



+ uX + v = 



(67) 



with 



u 



v 



2u 
27 




(68) 
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The nature of the roots is determined by the discriminant 



A = v 2 H u 3 . 

27 



(69) 



If A < we have case (i), if A = we have case (ii), and if A > we have 
case (iii). Thus, a 2 rlt is given by one of the two solutions of A = 0. These are 



It turns out that (a c J lt ) 2 ((«crit) 2 ) corresponds to the case where the imaginary 
parts of <7i and <T2 have different signs (the same sign, respectively). We con- 
clude that the physically interesting solution is (a^rit) 2 as gi ven by the first of 
Eqs. ( !70|) . For r = 1 (the two semicircles touch) it starts at unity and increases 
with r. 



10 Summary and Conclusions 

In this paper, we have defined and investigated two random-matrix models 
for the RPA equations. These models represent the generic forms of the RPA 
equations in cases of unitary or orthogonal symmetry. The RPA matrix is 
not Hermitean. Therefore, our random-matrix models do not belong to one 
of the ten classes listed in Ref. jX3] . Both random-matrix models depend on 
two parameters. The parameter a 2 gives the relative strength of the matrix 
C connecting states at positive and states at negative energies in relation to 
the strength of the matrix A, the matrix of interactions amongst the positive- 
energy states. For a 2 = 0, the average spectrum has the shape of two semicir- 
cles, one at positive and one at negative energies. The parameter r gives the 
distance of the centers of the two semicircles from the origin E = 0. 

We have focussed attention on the forms of the average spectra. We have 
not investigated the fluctuation properties of the positive- and the negative- 
energy branches of the spectra. This is because we are confident that both 
branches display the usual GUE or GOE level repulsion. Indeed, the total 
operator governing the positive-energy branch of the spectrum is given by 
A — 8 A where 6 A = C(E + r + A*)C* . The operator A — 8 A depends on 
energy. The addition of 5 A to A modifies the form of and shifts the average 
spectrum. But when we consider some energy interval comprising a finite 
number of levels in the limit iV — > oo, we may replace the energy argument E 
of 8 A by the value of E at the center of that interval and take 8 A as constant. 




(70) 
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The unitary or orthogonal invariance of the random matrix A — 5 A should 
then suffice to guarantee local spectral fluctuation properties of the GUE or 
GOE type. 

While level repulsion dominates the positive- and the negative-energy branches 
of the spectrum, there is level attraction between the two branches. That at- 
traction is due to the matrix C. The two semicircles describing the average 
spectrum for C = become deformed and move toward each other as the 
strength a 2 of C is increased. At some critical value a 2 = a^ rit , the spectra 
touch at E — 0. That value marks the point where the RPA equations become 
instable for most realizations of the random-matrix models. 

The mechanism that leads to instability is displayed by non-averaged RPA 
spectra. The symmetry of the RPA matrix implies that eigenvalues occur in 
pairs or in quartets: Real and purely imaginary eigenvalues occur pairwise 
with opposite signs; for fully complex eigenvalues with non-vanishing real and 
imaginary parts, the four eigenvalues lie symmetrically with respect to both 
the real and the imaginary axis in the complex E-plane. Level attraction 
leads to coalescence of pairs of real eigenvalues with opposite signs at E — 0. 
With increasing a 2 , the two eigenvalues leave the point E = and move in 
opposite directions along the imaginary axis. On that axis, pairs of eigenvalues 
may coalesce again. Because of the symmetry of the RPA equations, such 
coalescence must occur in parallel for the same value of a 2 on the positive 
and on the negative imaginary axis. A further increase of a 2 causes each pair 
to separate and to leave the imaginary axis in opposite directions along a 
straight line that runs parallel to the real E-axis. For very large values of a 2 , 
all eigenvalues come to lie on the imaginary E'-axis. The average spectrum 
has the shape of a semicircle centered at E — 0. 

We have studied these features both analytically and numerically. A gener- 
alized Pastur equation was derived for large matrix dimension (N — > oo). It 
reduces to a fourth-order polynomial for the spectral strength function <j\ of 
the positive-energy states. While for a 2 = ^sa\ is confined to positive en- 
ergies, So"i spreads to negative energies as a 2 is increased. We have shown 
that the level density obtained from the solution <j\ agrees very well with the 
result of numerical diagonalizations even of random RPA matrices of fairly 
low dimension. From the fourth-order equation for o\ we derived the critical 
value a 2 rit of the strength a 2 at which the two branches of the spectrum touch 
each other. That value is given by the simple equation a 2 rit = r 2 — 1 where r is 
the distance of the center of each semicircle from the origin E = measured in 
units of the radius of the semicircle. That result may be useful in applications: 
a 2 may be estimated by the ratio of two mean-square values measuring the 
strengths of the interaction matrix elements connecting particle-hole states 
at positive and negative energies (at positive energies only, respectively). The 
parameter r may be estimated as the ratio of the mean value of the positive 
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particle-hole energies and the root-mean-square value of the interaction ma- 
trix elements connecting states at positive energies. 
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Appendix 



We derive the eigenvalue decomposition (148]) of the trace of the Green's func- 
tion. In Chapter 1, Eq. (5.23) of Ref. [15] it is shown that the resolvent 
R{C) = {T — C) 1 °f an arbitrary matrix T in N dimensions with eigenvalues 
Xh, h = 1, 2, . . . , s possesses the decomposition 



m h -l 

R(0 = - E (c - a,)- 1 ^ + E (C - ^h)~ n ~ 1 D n 



h=l 



h 

n=l 



(71) 



Here ( is a complex variable and A^ ^ A/, for h ^ k. The matrices Ph are 
projectors which obey (Eq. (5.21) of Ref. [T3] ) 

P h P k = 5 hk P h ; E^ = 1 - (72) 

h 

The matrices D h obey (see Eqs. (5.24) and (5.50) of Ref. [13] ) 
P h D k = D k P h = 5 hk D k ; D h D k = (h ^ k) ; 

D™ h = where rrih = dimP^ . (73) 



In the case of Eq. (148]) . there are s = 2N distinct eigenvalues. This is true 
as long as the positive and negative energy spectra do not touch and is a 
consequence of level repulsion amongst the positive eigenvalues and amongst 
the negative eigenvalues. Therefore there are 2N different orthogonal projec- 
tors Ph which obey J^h^h — 1- This is possible only if dimP/j = 1 for all 
h = 1, 2, . . . , 2N. But m h = l for all h implies D h = for all h and Eq. (1711) 
reduces to 

27V 

^(CHE^-Cm. (74) 
h=l 

Since dimP^ = 1 we have TraceP^ = 1 for all h. Taking the trace of Eq. (1741) 
and identifying the variable ( with the energy, we obtain Eq. 
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